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Context. LS 5039 and LSI +61°303 are two binaries that have been detected in the TeV energy domain. These binaries are composed of a 
massive star and a compact object, possibly a young pulsar. The gamma-ray emission would be due to particle acceleration at the collision site 
between the relativistic pulsar wind and the stellar wind of the massive star. Part of the emission may also originate from inverse Compton 
scattering of stellar photons on the unshocked (free) pulsar wind. 

Aims. The purpose of this work is to constrain the bulk Lorentz factor of the pulsar wind and the shock geometry in the compact pulsar wind 
nebula scenario for LS 5039 and LSI +61°303 by computing the unshocked wind emission and comparing it to observations. 
Methods. Anisotropic inverse Compton losses equations are derived and applied to the free pulsar wind in binaries. The unshocked wind spectra 
seen by the observer are calculated taking into account the y — y absorption and the shock geometry. 

Results. A pulsar wind composed of monoenergetic pairs produces a typical sharp peak at an energy which depends on the bulk Lorentz factor 
and whose amplitude depends on the size of the emitting region. This emission from the free pulsar wind is found to be strong and difficult to 
avoid in LS 5039 and LSI +61 = 303. 

Conclusions. If the particles in the pulsar are monoenergetic then the observations constrain their energy to roughly 10-100 GeV. For more 
complex particle distributions, the free pulsar wind emission will be difficult to distinguish from the shocked pulsar wind emission. 



Key words, radiation mechanisms: non-thermal - 
X-rays: binaries 



stars: individual (LS 5039, LSI +61°303) - stars: pulsars: general - gamma rays: theory - 



o 

oo 
O 



1. Introduction 

Pulsars are fast rotating neutron stars that contain a large 
amount of rotational energy. A significant fraction of this en- 
ergy is carried away by an ultra-relativisti c wind of elec 



trons/positrons pairs and possibly ions (see iKirk et al.l 12007 



for a recent review). In the classical model of the Crab neb- 
' ula dRees & Gunnlll974l; iKennel & Coroniti|[l984l) . the pulsar 
wind is isotropic, radial and monoenergetic with a bulk Lorentz 
factor yo ~ 10^ far from the light cylinder where the wind is 
kinetic energy-dominated (cr <k 1). The cold relativistic wind 
expands freely until the ram pressure is balanced by the sur- 
rounding medium at the standoff distance Rs. In the termination 
shock region, the pairs are accelerated and their pitch angle to 
the magnetic field are randomized, producing an intense syn- 
chrotron source. Moreover, the inverse Compton scattering of 
the relativistic electrons on soft photons produces high energy 
(HE, GeV domain) and very high energy (VHE, TeV domain) 
gamma-rays. 

The shocked pulsar wind is thought to be responsible for 
most of the emitted radiation and gives clues about the proper- 
ties of this region. However, our knowledge of the unshocked 
pulsar wind region is limited and based on theoretical state- 



ments. If the magnetic field is frozen into the pair plasma 
as it is usually assumed, there is no synchrotron radiation 
from the unshocked wind. Nevertheless, nothing prevents in- 
verse Compton scattering of soft photons onto the cold ultra- 
relativistic pairs from occuring. The pulsar wind nebula (PWN) 
emission has two components: radiation from the shocked and 
the unshocked regions. 

Bogovalov & Aharonian ( 2000l) investigated the inverse 



Compton emission from the region upstream the termination 
shock of the Crab pulsar Comparisons between calculated and 
measured fluxes put limits on the parameters of the wind, 
in particular the si ze of the kinetic energy dominated region. 
iBall & KirkI (|2000|) investigated emission from an unshocked 
freely expanding wind with no termination shock in com- 
pact binaries. They computed spectra and light curves in the 
gamma-ray binary PSR B 1259-63, a system with a 48 ms pul- 
sar and a Be star in a highly eccentric orbit. The resulting 
gamma-ray emission is a line-like spectrum. 

In addition to PSR B1259-63, two other binaries have 
be en firmly confirme d as gamma-ray sources: LS 5039 



(lAharonian et al.ll2005l) and LSI H-6r303 dAlbert et al.ll2006l) . 
They are composed of a massive O or Be star and a com- 
pact object in an eccentric orbit. The presence of a young 
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pulsa r was detected only in PSR B 1259-63 jjohnston et al 
19921) . Radio pulses are detectable but vanish near perias- 
tron, probably due to free-free absorption and interaction with 
the Be disk wind. The compact PWN scenario is most prob- 
ably at work in this system and investigations w ere carried 
out to model high and very high energy radiation ( Kirk et al. 
1999 1: ISierpowska & Bednarekl 120051; 'Kh angulvan et al..,2007 : 
Sierpowska-Bartosik & Bednarek 2008). In LS 5039 and 
LSI +61°303 the nature of the compact object is still contro- 
versial but spectral and temporal similarities with PSR B1259- 
63 argue in fav or of the compact pulsar wind nebula scenario 
jPubusI 12006? ). The VHE radiation would therefore be pro- 
duced by the interaction between the pulsar wind and the stel- 
lar companion wind. The massive star provides a huge density 
of seed photons for inverse Compton scattering with the ultra- 
relativistic pairs from the pulsar wind. Because of the relative 
position of the compact object, the companion star and the ob- 
server, the Compton emission is modulated on the orbital pe- 
riod. The vicinity of a massive star is an opportunity to probe 
the pulsar wind at small scales. 

The component of the shocked pulsar wind was computed 
in lPubus et al. (2008) for LS 5039 and limits on the electron 
distribution, the pulsar luminosi ty and the magnetic field at the 
termination shock were derived. Sierp owska-Bartosik & TorresI 
(I2OO8.) calculated the VHE emission in LS 5039 as well, as- 
suming a power law injection spectrum for the pairs in the 
unshocked pulsar wind and pair cascading. In this paper, we 
investigate the anisotropic inverse Compton scattering of stel- 
lar photons on the unshocked pulsar wind within the compact 
PWN scenario for LS 5039 and LSI -i-61°303. Because of their 
tight orbits, the photon density is higher than in the Crab pulsar 
and PSR B 1259-63. A more intense gamma-ray signal from the 
unshocked pulsar wind is expected. The main purpose of this 
work is to constrain the bulk Lorentz factor yo of the pairs and 
the shock geometry. The next section presents the method and 
the main equations used in order to compute spectra in gamma- 
ray binaries. Section 3 describes and shows the expected spec- 
tra for LS 5039 and LSI +6r303 with diff'erent parameters. 
Section 4 discusses the spectral signature from the unshocked 
pulsar wind. 

2. Anisotropic Compton losses in y-binaries 

2.1. The cooling of pairs 

An electron of energy Ee = yenieC^ in a given soft photon field 
of density no ph cm"^ cools down through inverse Compton 
scattering (here the term 'electrons' refers indifferently to elec- 
trons and positrons). The power lost by t he electron is given by 
(|Joneslll965l:lBlumenthal & Gouldll970h 



the simple case of a monoenergetic and unidirectional beam of 
photons in the Thomson limit, the calculation of the Compton 
energy loss per electron is 



- ^ = o-Tcnoetiil -y6yUo)[(l -y6yUo)y^- 1] 



dt 



(2) 



where a-j is the Thomson cross section, jiQ - cos Oq and Oq 
the angle between the incoming photon and the direction of the 
electron motion. This calculation is done using the Compton 
kernel calculated by iFargion et al.l (Il997h . In the Thomson 
limit, the cooling of the electron follow a yl power law and 
has a strong angular dependance. In a more general way and 
for ye » 1 , the power lost per elec tron is calculated with the 
kernel derived in lDubus et al.l (l2008h Eq. (A. 6). 



2.2. Compton cooling of the free pulsar wind 

The pulsar is considered as a point-like source of monoener- 
getic and radially expanding wind of relativistic pairs e^/e". 
The pulsar wind momentum is assumed to be entirely caiiied 
away by the pairs. The companion star, with a typical luminos- 
ity of 10^** - 10^*^ ergs s"', provides seeds photons for inverse 
Compton scattering onto the radially expanding electrons from 
the pulsar The electrons see a highly anisotropic photon field. 
Inverse Compton efficiency has a strong dependence on do as 
seen in Eq. (|2]i. Depending of the relative position and direction 
motion of the electron with respect to the incoming photons di- 
rection, the cooling of the wind is anisotropic as well. Figure[T] 
sketches the geometry considered in the binary system to per- 
form calculations. 



Observer 



Pulsar 




dt j,_ 



eo)no——dei 
dtde\ 



(1) 



Companion star 



Fig. 1. Geometry of the binary system. Electrons of Lorentz 
factor ye are radially moving away at a distance r from the pul- 
sar and R from the companion star. The angle 41 quantifies the 
relative position between the pulsar, the companion star and the 
observer measures the angle between the electron direction 
of motion and the line joining the companion star center to the 
electron position through its motion to the observer 



where eo is the incoming soft photon energy, e\ the scattered 
photon energy and dN/dtdei is the Compton kernel. e± bound- 
aries are fixed by the relativistic kinematics of inverse Compton 
scattering. The cooling of the pairs e^/e" depends on the angu- 
lar distribution and spectrum of the incoming photon field. In 



For ultra-relativistic electrons , the radial dependence of the 
electron Lorentz factor jeif) for a given viewing angle is ob- 
tained by solving the first order differential equation Eq. ([U. 
iChernyakova & lUarionov ( 1999h found an analytical solution 
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in the Thomson Hmit and Ball & Kii5 ( 2000|) d erived a solu- 
tion in the general case using the IJonesI (Il965h results for a 
point-like and monoenergetic star with » 1 . In this approxi- 
mation, the density of photons is L*/(47rc/?^eo) ph cm"^, where 
is the star luminosity and eo - 2.7A;r* the average energy 
photon from the star The differential equation is then 



dr 



1 



rrieC^ AncE} 



dN 
dtde\ 



dei 



(3) 



where - d^ + - Irdcosiff. Calculations beyond the mo- 
noenergetic and point-like star approximation require two extra 
integrations, one over the star spectrum and the other onto the 
angular distribution of the incoming photons due to the finite 
size of the star The complete differential equation is then given 
by 



dr nieC^ Iff ^ ' 



dN 

■ eo) »o , , deideodD.0- 
atae\ 



(4) 



For a blackbody of temperature T^, and a spherical star of ra- 
dius Rj,, the incom ing photon density hq is given by Eq. (13) in 



Dubus et alj (120081) . It is more convenient to compute the calcu- 
lation of the Lorentz factor as a function of rather than r (see 
Fig. [ill. These two variables are related through the relation 

/ tan (A \ 

r^c/cosi/r 1 , re[0, -Hoo], i/^^ e [i/r, tt] . (5) 

\ tani/r,./ 

Figure|2]presents the numerical computed output solution yiij/,) 
applied to LS 5039 with an inclination of / = 60° for a neutron 
star where the viewing angle varies between 7r/2 - / = 30° and 
njl + i = 150°. Here, the wind is assumed to have an injection 
Lorentz factor y(i^, (0)) = jo - 10^ and to continue unimpeded 
to infinity (i.e. it is not contained by the stellar wind). 

For small viewing angles ifr, the cooling of the wind is 
very efficient because the collision electron/photon is almost 
head-on and the electrons are moving in the direction of the 
star where the photon density increases. For viewing angles 
t// ^ ;r/2, the cooling of the pairs is limited. In all cases, most 
of the cooling occurs at ifr^ ~ lA- For i^,. n/2, the electron 
is moving away from the star and the scattering angle become 
small leading to a decrease in the wind energy loss. A com- 
parison of Compton cooling between the point-like and finite 
size star is shown in Fig.|2] The effects of the finite size of the 
star are significant in two cases. The impact of the finite size of 
the star is important if the observer is within the cone defined 

I 1! — — ! 

by the star and the electron at apex (see Dubus et al. 2008 for 
more details). For viewing angles iff ^ arcsin(/?*/(i), the cool- 
ing is less efficient whereas for tfr ^ n - arcsin(/?* /c/) it is more 
efficient as it can be seen in the two extreme value of i// in figure 
|2l The other situation occurs when the electrons travel close to 
the companion star surface, for i// ^ jt/2 and tf^r ^ n/2. In that 
case the angular distribution of the stellar photons is broad and 
close head-on scatterings are possible, leading to more efficient 
cooling compared with a point-like star. Nevertheless, these ef- 
fects remain small for LS 5039 and LSI +61°303 and will be 
neglected in the following spectral calculations. 
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Fig. 2. Compton cooling of a monoenergetic, free pulsar wind 
with 70 = 10^ d = 2R^ (T^ = 39 000 K, R^ = 9.3 Rq). The 
different curves show the dependence with the viewing angle ifr 
on the cooling, if/ varies between 30° (bottom) and 150° (top) if 
/ = 60°. Each curve shows the evolution of the Lorentz factor 
y with if/r as the electron moves along the line of sight, if/r is 
related to r by Eq. (|5]l so that if/,- = if/ at r = Q and if/r = tt for 
r - +oa. The calculation was carried out for a blackbody point 
like star (solid line) and taking into account the finite size of 
the star (dashed line). 

2.3. Unshocked pulsar wind spectra 

The number of scattered photons per unit of time, energy and 
solid angle depends on three contributions: the density of the 
incoming photons, the density of target electrons and the num- 
ber of scattered photons per electron. The pulsar wind of lumi- 
nosity Lp is assumed isotropic and monoenergetic, composed 
only of pairs and with a negligible magnetic energy density 
(cr <K 1). The electrons density (e" cm erg ') is then propor- 
tional to Ifr^ if pair production is neglected. Here, the interest- 
ing quantity for spectral calculations is the number of electrons 
per unit of soUd angle, e nergy and length, w hich is time the 
electrons density so that jBall & Kirkll2000l) 



dNe 



dQ.edydr Anc/SoyomeC 



■S(y-ye(r)), 



(6) 



with 6 the Dirac distribution. In deep Klein-Nishina regime, 
spectral broadening is expected because the continuous energy 
loss prescription fails (AEe ~ Eg). The complete kinetic equa- 
tion must be u sed in order to describe acc urately the electrons 
dynamics (see lBlumenthal & Gouldlfl970l Eq. (5.7)). However, 
the 6 approximation used here is reasonably good (Zdzia rskl 
1989b . The case of an anisotropic pulsar wind is discussed in 
§4. In the following, the pulsar wind will therefore be assumed 
to follow Eq. (|6]l. Heating of the pulsar wind by the radiative 
drag is neglected (Ball & Kirk 200Q). 

In order to compute spectra, the emitted photons are sup- 
posed to be entirely scattered in the direction of the electron 
motion. Because of the ultra-relativistic motion of the elec- 
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Fig. 3. Computed inverse Compton spectrum from the unshocked pulsar wind in LS 5039 and its dependence with the emitting 
region size Rg. The pulsar wind has yo - 10^, Lp = 10^^ erg s"' and the star is a point-like blackbody. Spectra are calculated at 
the superior (left) and inferior (right) conjunctions for different standoff distances Rs - 10'" (bottom), 3 10'", 10", 3 10" cm 
and +00 (dashed line). 



trons, most of the emission is within a cone of aperture angle of hydrodynamical limit. They obtained the geometry for the rela- 



the order of 1 /y^, <K 1 . In this classical approximation, the spec- 
trum seen by the observer is the superposition of the contribu- 
tions from the electrons along the line of sight pulsar-observer 
in the solid angle dQ.e. The spectrum seen by the observer is 
obtained with the following formula 



dNtc 



dtde\ dQ.1 



-M 



no 



dN 



dNe 



dtde\ dQ.edydr 



dyd£odrdQ.Q (7) 



where Tyy takes into account the absorption of gamma-rays due 
to pair production with soft photons from the companion star 
and is calculated following , Dubua (.2006a.) . 



2.4. The compact PWN geometry 

The collision of the relativistic wind from the pulsar and the 
non-relativistic wind from the massive star produces two termi- 
nation shock regions separated by a contact discontinuity (see 
Fig.lDi. The geometry of the shock fronts are governed by the 
ratio o f the fl ux wind momentum quantified by 77 and defined 
as re.g. lStevens et aL.1992;.Eichler & Usov. 19931) 



(8) 



where M„ is the mass loss rate and Voo the stellar wind speed 
of the O/Be star. For two spherical winds, the standoff" distance 
point Rs depends on 77 and on the orbital separation d 



1 + 



(9) 



Bogovalov et alj (l200a) have investigated the collision between 
the pulsar wind and the stellar wind in the binary PSR B1259- 
63, with a relativistic code and an isotropic pulsar wind in the 



tivistic and nonrelativistic shock fronts and the contact discon- 
tinuity. They find that the collision between the two winds pro- 
duces an unclosed pulsar wind termination shock (in the back- 
ward facing direction) for 77 > 1.25 10"^. 



Observer 




Companion star 



Fig. 4. Shock geometry considered for the wind collision. For 
rj > 1.25 10"^, the pulsar wind region remains open with 
an asymptotic half-opening angle a. The dark region is the 
shocked relativistic pulsar wind and the light region is the 
shocked non-relativistic stellar wind, separated by a contact 
discontinuity (dot-dashed line). The size of the emitting zone 
seen by the observer depends on the viewing angle ^. 



The size of the emitting region depends on the shock ge- 
ometry and the viewing angle, whi ch can therefo r e have a ma- 
jor impact on the emitted spectra. iBall & Doddl (1200 lb com- 
puted spectra from the unshocked pulsar wind in PSR B1259- 
63 for an hyperbolic shock front terminated close to the pulsar. 
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They found a decrease in the spectra fluxes and a decrease in 
the hght curve asymmetry and flux parti cularly near perias tron 
compared with the spectra computed bv lBall & KirkI ( 120001) . 

Figure |3] presents computed spectra, ignoring 77 absorp- 
tion at this stage, applied to LS 5039 at the superior and infe- 
rior conjunctions for diff'erent standoff distances and a pul- 
sar wind with yo - 10''. At the superior conjunction where 
t// = 30°, the Compton cooling of the wind is eflicient. The 
broadness in energy of the radiated spectra is related to the size 
of the unshocked pulsar wind region. For small standoff dis- 
tances Rs d, spectra are truncated and sharp because the 
termination shock region is very close to the pulsar, so that the 
pairs do not have time to radiate before reaching the shock. 
For Rs ^ d, the free pulsar wind region is extended and emis- 
sion from cooled electrons starts contributing to the low energy 
tail in the scattered spectrum. The amplitude of the spectrum 
reaches a maximum when the injected particles can cool effi- 
ciently before reaching the shock. The spectral luminosity is 
then set by the injected power and is not affected anymore by 
the size of the emitting zone. At the inferior conjunction where 
t// - 150°, the cooling is less efficient and most of the emis- 
sion occurs close to the pulsar where the photon density and Oq 
are greater, regardless of the size of the emitting region. The 
radiated flux then depends linearly on R,. A complete inves- 
tigation is presented in the next section where absorption and 
spectra along the orbit are computed and applied to LS 5039 
and LSI +61°303, ignoring pair cascading. 

3. Spectral signature of a monoenergetic pulsar 
wind in LS 5039 and LSI +6r303 

In the following sections, the emission expected by the un- 
shocked pulsar wind in LS 5039 and LSI -h61°303 is compared 
with measured fluxes. Because spectra depends on the shock 
geometry and the injection Lorentz factor, spectra are calcu- 
lated for various values of the two free parameters rj and 70. 

3.1. LS5039 

The companion star and the pulsar winds are assumed isotropic 
and purely radial. The orbital par ameters are those m easured by 
Casares et al (l2005bh as used in lDubus et al.l (l2008h . 

Figure |5] gives the total power radiated by the electrons 
in the unshocked pulsar wind as a function of 70 at peri- 
astron. Here, the shock front is assumed spherical of radius 
Rs- A maximum of efficiency is observed at about 70 ~ 10^ 
which corresponds to the transition between the Thomson and 
Klein-Nishina regimes where the Compton timescale is short- 
est (Dubus 2006 b). The fraction of the pulsar wind power ra- 
diated at periastron depends strongly on Rs- It is about 20% 
for 77 = 10"^ and can reach 70% for rj = 0.1. Hence, most of 
the spindown energy can be radiated directly by the unshocked 
pulsar wind. 

Figure |6] presents computed spectra averaged along the or- 
bit for different shock geometry and Lorentz factor with a pul- 



sar spindown luminosity of Lp 



10^^ erg s-\ 



The relativistic 
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Fig. 5. Total radiated power by the unshocked pulsar wind Prad 
in LS 5039 as a function of 70. Prad is computed at periastron 
for 7/ = 10"^ (solid line), 2 10"^ (dotted Hne), 10"' (dashed 
line) and with no termination shock (dotted-dashed line). 



angle a is taken from Eq. (27) in lBogovalov et al.l (120081) . both 
parameters depending only on 77. Figure |4] sketches the shock 
morphology for 1.25 10"^ < 77 < 1 and presents the different 
shock fronts expected. The twist due to the orbital motion is 
ignored since most of the emission occurs in the vicinity of the 
pulsar. The size of the emitting zone seen by the observer is 
thus set for any given viewing angle if). Note that it is always 
greater than R,. The remaining free parameter 70 is chosen in- 
dependently between lO'* and 10^. 

Computed spectra predict the presence of a narrow peak 
in the spectral energy distribution due to the presence of the 
free pulsar wind. The luminosity of this narrow peak can be 
comparable to or greater than the measured fluxes by EG RET 



wind parameters to ~ 10 ^ Mq yr ' and Voo ~ 2400 km s ' 

lO-'^ erg s"' then gives 



shock front is described by an hyperbolic equation. The hy- 
perbola apex is set by Eq. (|9]l and the asymptotic half-opening 



and HESS dHartman et al.1 119991: 1 Aharonian et al .1120061) . For 
J] - 10"^, the pulsar wind termination shock is closed and 
the unshocked wind emission zone is small. For 77 = 0.02 and 
?7 = 0.1 the line spectra are well above both the limits imposed 
by the HESS observations. The extreme case with no termina- 
tion shock shows little differences with the case where 77 = 0.1. 
Spectroscopic observations of LS 5039 constrains the O star 

J ■■ ^1 k 

(iMcSwain et al.ll2004l) . Assuming Lp 

77 ~ 2 10' ^ (top right panel of Fig. |6]l or Rs ~ 2 lO" cm as 
m (lDubusl[2006b ). In this case, almost half of the pulsar wind 
energy is lost to inverse Compton scattering before the shock is 
reached (Fig.|5]l. This is an upper limit since the reduced pulsar 
wind luminosity would bring the shock location closer to the 
pulsar than estimated from Eq. (|9]l. HESS observations already 
rule out a monoenergetic pulsar wind with 70 = 10^ or 10^ and 
Lp = lO^*" erg s"' as this would produce a l arge component eas - 
ily seen at all orbital phases (see Fig. 5 in lDubus et al.ll2008h . 
The EGRET observations probably also already rule out values 
of7< 10^. 




Fig. 6. Spectral signature from the unshocked pulsar wind expected in LS 5039 and dependence with yo and 77. Spectra are 
averaged on the orbital phases corresponding to the HESS 'high state' (solid line, 0.45 < (p < 0.9, with <p = at periastron) and 
'low state' (dashed line, (p < 0.45 or (p > 0.9). The spectra are compared with EGRET (dark bowtie) and HESS (light bowties) 
observations, adopting a distance of 2.5 kpc. In the top left panel, 77 = 10"^ the shock is closed and the unshocked pulsar wind is 
assumed spherical. For tj - 2 10"^ (top right panel) and 77 = 0.1 (bottom left panel) the shock is open with half-opening angles 
a ~ 2° and a ~ 30° respectively. The bottom right panel shows the extreme case with no termination shock. 



3.2. LSI +61° 303 

In this system the stellar wind from the companion star is as- 
sumed to be composed of a slow dense equatorial disk and 
a fa st isotropic polar win d. The stellar wind may be clumpy 
and lZdziarski et al.l (l2008b have proposed a model of the high- 
energy emission from LSI +61°303 that entails a mix between 
the stellar and the pulsar wind. The o rbital parameters are those 
measured by ICasares et al (l2005ah (new orbital parameters 
were recently measured by lGrundstrom et al ] l2007h . 



Computed spectra applied to LSI -i-61°303 and averaged 
over the orbit to compare with EGRET and M AGIC luminosi- 
ties jHartman et al.lll999l: IXlbert et al.ll2006h are presented in 
Fig. |8] New data were recently reported by the MAGIC collab- 
oration (1 Albert et al., 2008) . They confirmed the measurements 
of the first observational compaign and found a periodicity in 
the gamma-ray flux close to the orbital period. The pulsar spin- 
down luminosity is set to Lp - 10^^ erg s ' and the injected 



Lorentz factor to 10^ 10^ 10^ and lO'' as for LS 5039. There is 
more uncertainty in rj because of the complexity of the stellar 
wind. The polar outflow is usuall y modelled with M „ = 10 
Mq yr ' and Voo = 2000 km s ' d Waters et al.ll 19881) leading to 
7] ~ 0.2 - 0.3. Concerning the slow dense equatorial disk, the 
mass flux is typically one hundred times greater than the polar 
wind and the terminal velocity is a few hundred km i"' giving 
T] compatible with ~ lO"-' - 10"^. 

The overall behaviour is similar to LS 5039. The spectral 
luminosities and the total power radiated by the unshocked pul- 
sar wind (Fig.© are lower in LSI -H6r303 than LS 5039 be- 
cause the compact object is more distant to its companion star 
and the latter has a lower luminosity, leading to a decrease in 
the density of seed photons for inverse Compton scattering. If 
7] - lO"-', no constrains on yo can be formulated as the spec- 
trum is always below the observational limits. For larger val- 
ues of 77, the very high energy observations constrain yo to be- 
low 10^, assuming the pulsar wind is monoenergetic. Spectra 
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Fig. 7. Total radiated power by the unshocked pulsar wind P^ad 
in LSI +61°303 as a function of yo. P,.,,^ is computed at perias- 
tron for 77 = 10"^ (solid line), 2 10"^ (dotted line), 0.53 (dashed 
line) and with no termination shock (dotted-dashed line). 
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and 



were computed for 77 = 0.5 3 with M„ = 10 
Voo = 1000 km s ' as used bv lRomero et alj(l2007h . In this case. 



the spectra are close to the freely propagating pulsar wind. The 
EGRET luminosity is slightly overestimated for yo < 10^ when 
77 = 0.53. 

4. Discussion 

The proximity of the massive star in LS 5039 and LSI +61°303 
provides an opportunity to directly probe the distribution of 
particles in the highly relativistic pulsar wind. The calcula- 
tions show the inverse Compton emission from the unshocked 
wind should be a significant contributor to the observed spec- 
trum. For a monoenergetic and isotropic pulsar wind the emis- 
sion remains line-like, with some broadening due to cooling, 
as had been found previously for the Crab and PSR B 1259-63 
dBogovalov & AharonianlboOOl: iBall & KirklboOOh . However, 
here, such line emission can pretty much be excluded by the 
available very high energy observations of HESS or MAGIC, 
and (to a lesser extent) by the EGRET observations that show 
power-law spectra at lower flux levels. 

4.1. Is the pulsar wind power overestimated? 

Reducing the pulsar power (or, equivalently, increasing the dis- 
tance to the object) would diminish the predicted unshocked 
wind emission relative to the observed emission. This is not 
viable as this would also reduce the level of the shocked pul- 
sar wind emission. Similarly, the energy carried by the parti- 
cles may represent only a small fraction of the wind energy. 
At distances of order of the pulsar light cylinder the energy 
is mostly electromagnetic. Evidence that this energy is con- 
verted to the kinetic energy of the particles comes from pleri- 
ons, which probe distances of order 0. 1 pc from the pulsar. It 



is therefore conceivable that this conversion is not complete at 
the distances under consideration here (0.01-0.1 AU). In this 
case the emission from the particle component would be re- 
duced. However, the shocked emission would also be reduced 
as high cr shocks divert l ittle of the energy into the particles 
dKennel & Coronitilll984 . ] -furthermore, the high energy parti- 
cles would preferentially emit synchrotron rather than inverse 
Compton due to the higher magnetic field. Hence, this possi- 
bility also seems unlikely. 

Alternatively, the unshocked wind emission could be 
weaker compared to the shocked wind emission if the termi- 
nation shock was closer to the pulsar, i.e. if one had a low 
77. In LS 5039 the unshocked wind emission is strong even 
with 77 - 0.001, which already implies a stronger stellar wind 
than optical observations seem to warrant. Furthermore, the 
value of the magnetic field would be high if the termination 
shock was close to the pulsar and this inhibits the formation 
of very high energy gamma-rays as the high energy electrons 
would then preferentially lose energy to synchrotron radiation 
(Dubus 2006b). Hence, it does not seem viable either to invoke 
a smaller zone for the free wind. 

The conclusion is that the strong emission from the pul- 
sar wind found in the previous section is robust against general 
changes in the parameters used. The following subsections ex- 
amine how this emission can be made consistent with the ob- 
servations. 



4.2. Constrains on the pulsar wind Lorentz factor 

The high level of unshocked emission is compatible with the 
observations only if it occurs around 10 GeV or above 10 TeV, 
i.e. outside of the ranges probed by EGRET and the current 
generation of Cherenkov telescopes. This poses stringent con- 
straints on the energy of the particles in the pulsar wind. The 
Lorentz factor of the pulsar wind would be constrained to a 
few 10^ or to more than 10^. The 1-100 GeV energy range 
will be partly probed by GLAST and HESS-2, and CTA in the 
more distant future. For instance, unshocked wind emission in 
LS 5039 would appear in the GLAST data as a spectral harden- 
ing at the highest energies. Nevertheless, that the free wind is 
emitting in the least accessible spectral region may appear too 
fortuitous for comfort. 



4.3. Anisotropic pulsar wind 



The assumptions on the pulsar wind may be inaccurate. 
Pulsa r winds are thought to be anisotropic (Begel man & Lil 



1992h . iBogovalov & Khangoulvan (2002) interpreted the jet- 



torus structure revealed by X-ray Chandra observations of the 
Crab nebula, as a latitude dependence of the Lorentz factor 
j{0) - ji + y„, sin^ 6 where y; is small (say 10^) and y„, is high 
(say 10^). This hypothesis was corroborated by computational 
calculations in Komissar ov & Lvubarskvl (1200 40 where the syn- 
chrotron jet-torus was obtained. Here, the pulsar orientation 
to the observer is fixed (unless it precesses) so that the initial 
Lorentz factor of the pulsar wind along the line of sight would 
remain the same along the orbit. However, assuming the par- 




Fig. 8. Spectral signature from the unshocked pulsar wind expected in LSI +61° 303 and dependence with jq and rj. Spectra are 
averaged between phase 0.4 < <p < 0.7 (solid line, at periastron <p = 0.23) and the complementary phases < 0.4 or > 0.7 
(dashed line). Luminosities are compared with the EGRET (dark bowtie) and MAGIC (light bowtie) observations, adopting a 
distance of 2.3 kpc. In the top left panel, rj = 10"^ the shock is closed and the unshocked pulsar wind is assumed spherical. For 
T] - 2 10"^ (top right panel) and 77 = 0.53 (bottom left panel) the shock is open with half-opening angles 0-^2° and a ~ 60° 
respectively. The bottom right panel shows the extreme case with no termination shock. 



tide flux in the pulsar wind remains isotropic, the unshocked 
wind emission will appear at a lower energy and at a lower 
flux if the pulsar is seen more pole-on. The peak energy of the 
line-like spectral feature directly depends on y(0). Its intensity 
will also decrease in proportion as the pulsar power matches 
the latitude change in y to keep the particle flux isotropic (see 
Eq. 6). 

The shocked wind emission is set by the mean power and 
Lorentz factor of the wind and is insensitive to orientation. 
However, a more pole-on orientation will lower the contribu- 
tion from the unshocked component. For instance, if y{6) = 
lO'* + lO'^sin^e and 6 - IT.5 then the effective y along the 
line-of-sight will be 10^ and the observed luminosity of the un- 
shocked emission will be lowered by a factor 10 compared to 
the mean pulsar power (Eq. 6). The probability to have an ori- 
entation corresponding to a value of y(0) of 0. 1 y,,, or less is 
about i, assuming a uniform distribution of orientations. This 
would be enough to push the line emission to lower energies 



and to lower fluxes by a factor 10 or more, thereby relaxing the 
constraints on the mean Lorentz factor of the wind. Although 
this is not improbable, it would again require some fortuitous 
coincidence for the pulsars in both LS 5039 and LSI -i-61°303 
to be seen close enough to the pole that their free wind emission 
is not detected. 

4.4. The energy distribution of tine pairs 

The assumption of a monoenergetic wind may be incorrect, if 
only because the particles in the pulsar wind are bathed by a 
strong external photon field even as they accelerate and that 
this may lead to a significantly different distribution. Fig. |9] 
shows the emission from a pulsar wind where the particles have 
been assumed to have a power-law distribution with an index 
of -2 between y of lO^* and 10^. Obviously, a power-law distri- 
bution of pairs erases the line-like spectral feature. The emis- 
sion properties are essentially identical to the emission from 
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the shocked region with a harder and fainter intrinsic Compton 
spectrum when the pulsar is seen in front of the star compared 
to when it is behind. The emis sion from th e shocked region 
is also shown, calculated as in iDubus et al. I (12008,) . The par- 
ticle injection spectrum is the same in both regions. The par- 
ticles are assumed to stay close to the pulsar and to escape 
from the shocked region after a time fgsc = Rs/ic/3) (top) and 
10/?s/(c/3) (bottom). Longer fesc do not change the distribution 
any further. The longer escape timescale enables a harder parti- 
cle distribution to emerge at high energies (where the radiative 
timescale is comparable to Rs/c, see Fig. 2 in Dubus 2006b). 
With fesc - Rs/ic/3), the shocked spectra is very close to the 
unshocked spectra. Generally, calculations show the spectra 
from the shocked and unshocked regions may be indistiguish- 
able when the injected particles are taken to be the same in 
both regions. The only possible difference is that the longer res- 
idence time of particles in the shocked region allows for harder 
spectra. 



4.5. Dominant emission from the unshocked wind 

Emission from the unshocked wind could be the dominant con- 
tributor to the spectral energy distribution. In this case, the 
observations give the particle distribu tion in the pulsar wind. 
Sierpowska-Bartosik & TorresI (l2008h have considered such a 
scenario for LS 5039 and use a total energy in leptons of 
about lO"*^ erg s ' and a power-law distribution with an index 
around -2, both of which are adjusted to the observations and 
vary with orbital phase. The total pulsar power is much larger. 



10 ,nll mis inl3 



10 erg s , in order to have a big enough free wind emis- 
sion zone. Most of the pulsar energy is then carried by nuclei. 
Such a large luminosity would make the pulsar very young, 
comparable to the Crab pulsar, implying a high birth rate for 
such systems. Fig.fTOlshows the expected emission from a pul- 
sar wind propagating to infinity and with a particle power law 
index of -L5 from y - 10^ to 10^ chosen to adjust the 'high 
state' of LS 5039. The injected power is 4 10^^ erg s"'. The 
injected spectrum gives a good fit of the 'hard' state. However, 
the 'low' state dominates the complete very high energy contri- 
bution (^ 1 TeV) due to the extended emitting region. Particles 
have enough time to radiate very high energy gammay-ray far 
away, where they are less affected by y - y absorption. 

A possible drawback of this scenario is that the synchrotron 
and inverse Compton emission are not tied by the shock condi- 
tions and that the total energy in leptons is low so that it is not 
clear how the radio. X-ray and gamma-ray observations below 
a few GeV would arise. It is also unclear how this can lead to a 
coUimated radio outflow as seen in LSI +61°303. A possibility 
is emission from secondary pairs created in the s t ellar w ind by 
cascading as suggested by Bosch-Ramon et al. (12008V More 
work is necessary to understand these different contributions 
and the signatures that may enable to distinguish them. 

Another potential drawback of this scenario is that it does 
not explain why the very high energy gamma-ray flux is ob- 
served to peak close to apastron in LSI -i-61°303. If the inverse 
Compton scattering in the pulsar wind is responsible then the 
maximum should be around periastron, especially as gamma- 
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Fig. 9. Comparison between emission from the shocked and 
unshocked regions in LS 5039, taking the same particle injec- 
tion for both regions. The distribution is a power-law of index 
-2 between 7=10"' and 10*^ with total power lO"'^ erg s"'. 
Spectra are averaged to correspond to the HESS 'high state' 
(solid lines) and 'low state' (dashed lines) as in Fig. |6] The 
geometry is a sphere of radius 2 10" cm. Grey lines show 
emission from the unshocked emission and dark lines show the 
emission from the shocked region. The orbital averaged non- 
absorbed spectrum from the unshocked pulsar wind is shown 
in grey dotted line. Particles escape from the shocked region on 
a timescale fesc = Rs/(c/3) (top) or lORs/(c/3) (bottom). The 
unshocked emission is the same in both panels. 



gamma attenuation is very limited in LSI -i-61°303. On the 
other hand, if the emission arises from the shocked pulsar wind 
then synchrotron losses explain the lack of very high energy 
gamma-rays at periastron: the pulsar probes the dense equa- 
torial wind from the Be star and the shock forms at a small 
distance, leading to a high magnetic field and a cutoff in the 
high-energy spectrum (see §6.2.2 in iDubus..2006b,) . 
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Fig. 10. Fit to the observations of LS 5039, assuming the domi- 
nant contribution comes from the unshocked region of the pul- 
sar wind. The 'high state' of HESS corresponds to the solid 
line, the 'low state' to the dashed line and the dotted line to the 
orbital averaged non-absorbed spectrum. The injected particles 
had a power law distribution with an index of - 1 .5 between 
y = lO-* and 10**. The injected power is 4 10^^ erg s~'. 



5. Conclusion 

Gamma-ray binaries are of particular interest in the study of 
pulsar wind nebula at very small scales. The massive star radi- 
ates a large amount of soft seeds photons for inverse Compton 
scattering on relativistic electrons from the pulsar. One expects 
two contributions in the gamma-ray spectral energy distribu- 
tion: one from the shocked pulsar wind and another from the 
unshocked pulsar wind. The spectral signature from the un- 
shocked region is strong and depends on the shock geometry 
and the initial energy of the pairs in the pulsar wind. A signif- 
icant fraction of the pulsar wind power can be lost to inverse 
Compton scattering before a sh ock forms with the stellar wind 
(iSierpowska & Bednarekll2005h . The shock location will thus 
be slightly closer in to the neutron star than calculated without 
taking into account the losses in the pulsar wind, assuming the 
wind is composed only of e^e" pairs. 

Significant emission from the free pulsar wind seems un- 
avoidable. Inverse Compton losses in the free wind may be re- 
duced if the shock occurs very close to the neutron star. This 
is unlikely as it would require a very strong stellar wind or a 
pulsar wind power that would be too low to produce the ob- 
served emission. If the pulsar wind is anisotropic then the ori- 
entation of the pulsar with respect to the observer can make the 
unshocked emission less conspicuous. This comes at the price 
of a peculiar orientation. If the pulsar wind is monoenergetic, 
then the line-like expected spectrum exceeds the observed very 
high energy power-laws for all geometries unless the pair en- 
ergy is around 10 GeV or above 10 TeV. This pushes the direct 
emission from the wind in ranges where it may be constrained 
by future GLAST, HESS-2 or CTA measurements. The absence 



of any line emission will show that the assumption of a Crab- 
like monoenergetic, low cr pulsar wind was simplistic. If the 
pairs in the pulsar wind have a power-law distribution, then 
the unshocked emissi on is essentially indistinguishable fr om 
the shocked emission ( Sierpowska-Bartosik & Torresll2008 ). A 
promising alternative is the striped wind model in which the 
wind remains highly magnetised up to the termination shock, 
where the alt ernating field cou ld be dissipated and accelerate 
particles (see iKirk et al J 120071 and references therein). Future 
theoretical studies on the generation of pulsar relativistic winds 
in the context of a strong source of photons may be able to yield 
the particle distribution to expect and lead to more accurate pre- 
dictions. 
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